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Abstract 

We compare recent lattice studies of QCD thermodynamics at non-zero quark 
chemical potential with the thermodynamics of a hadron resonance gas. We argue 
that for T < T c the equation of state derived from Monte-Carlo simulations of 
two flavour QCD at non-zero chemical potential can be well described by a hadron 
resonance gas when using the same set of approximations as used in current lattice 
calculations. We estimate the importance of truncation errors arising from the use 
of a Taylor expansion in terms of the quark chemical potential and examine the 
influence of unphysically large quark masses on the equation of state and the 
critical conditions for deconfinement. 
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1 Introduction 



While the thermodynamics of strongly interacting matter at vanishing baryon 
number density or chemical potential has been studied pQ in lattice calcu- 
lations for quite some time, the first investigations of the equation of state 
at non-vanishing quark chemical potential (fi q ) have started only recently 
[21 El El- These studies of bulk thermodynamics have been performed with 
different lattice actions and also have used different methods (exact ma- 
trix inversion [2] or Taylor expansion [3E]) to extend previous calculations 
performed at n q — into the region of \i q > 0. Nonetheless they led to 
qualitatively and even quantitatively similar results. 

The basic pattern found for the temperature dependence of, e.g. the 
pressure at finite density, follows closely that already seen at \i q = 0; the 
pressure changes rapidly in a narrow temperature interval and comes close 
to the Stefan-Boltzmann value of an ideal gas of quarks and gluons at about 
twice the transition temperature. Consequently, the density dependence of 
the equation of state in the high temperature plasma phase has successfully 
been compared [S] with quasi-particle models that were also used at \i q = 0. 

In this work we concentrate on a discussion of the thermodynamics of the 
hadronic phase of QCD in the regime of low baryon number density (fi q /T<l) 
but high temperature (T ~ T c (fi q = 0)). In a recent paper we have shown jB] 
that the partition function of a hadron resonance gas yields quite a satisfac- 
tory description of lattice results on bulk thermodynamic observables in the 
low temperature, hadronic phase of QCD at n q = 0. We will extend here our 
previous study to finite chemical potential and compare the predictions of 
the resonance gas model calculations with lattice results obtained in 2-flavour 
QCD using a Taylor expansion for small [i q fT 0]. The reference system for 
these calculations is a previous analysis [7j of the temperature dependence 
of the pressure in 2-flavour QCD performed at fi q = 0. Unlike the approach 
based on an exact inversion of the fermion determinant [2] the Taylor ex- 
pansion, obviously, has the disadvantage of being approximate. There is, 
however, good reason to expect, that at least at high temperature, the con- 
tribution of terms that are beyond (9((/x g /T) 4 ) order is small. Moreover, as 
will become clear from our discussion below, it turns out that the expansion 
coefficients themselves provide useful information on the relevant degrees of 
freedom controlling the density dependence of thermodynamic quantities. 
We argue that baryons and their resonances are these relevant degrees of 
freedom that govern thermodynamics in the confined phase at finite density. 
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We show that for T <T C the equation of state at non-zero chemical potential 
which has been obtained in lattice calculations can be well described by a 
baryonic resonance gas when using the same set of approximations as used 
in current lattice studies. We examine the importance of truncation effects 
in the Taylor expansion and discuss the influence of unphysically large quark 
mass values on thermodynamic observables and the critical conditions for 
deconfmement. 



2 Finite density QCD and Taylor Expansion 

The basic quantity that describes thermodynamics at non-vanishing chemical 
potential is the pressure. In the grand-canonical ensemble it is obtained 1 
from the partition function as 

T 

p(T, Hq) = Jim — In Z(T, fj, q , V) . (1) 
For small values of n q /T the pressure may be expanded in a power series, 

. (2) 

In recent lattice calculations this series has been analyzed up to C((/i 9 /T) 4 ) 
and in addition to the density dependent change of the pressure (Ap) quanti- 
ties like the quark number density (n q ) and quark number susceptibility (x q ) 
have been calculated. The latter are obtained from Eq.|2]through appropriate 
derivatives with respect to the quark chemical potential, 

Ap = p{T, N ) -p(T,0) (fh\ a + CA(T) . 



n q 
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T 3 dfi q \T J " v ' \ T 

Xq 2 p(T,H q ) rrr . frr . (H 



T 2 dill y ' y ' \ T 



2 c 2 (T) + 12 c 4 (T) ^ • (3) 



1 Although the volume dependence of thermodynamic quantities calculated on the lat- 
tice requires a careful analysis and has not yet been performed for most thermodynamic 
studies with dynamical quarks, we will for simplicity of notation suppress in the following 
any volume dependence in our formulas. 
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In the asymptotically high temperature limit this expansion terminates, in 
fact, at the order given in Eq. El The expansion coefficients are then given 
by €2(00) = Uf/2 and 04(00) / '02(00) = l/2n 2 respectively 2 . In this limit the 
ratio C4/C2 is small and the leading order term, consequently, dominates in 
the Taylor expansion for a wide range of values for fi q /T. The lattice results 
for the expansion coefficients, obtained in 2-flavour QCD are shown in Fig. ^ 
It can be seen, that at T ~ 1.5 T the numerical results for &i(T) still deviate 
by about 20% from the ideal gas value while the ratio C4/C2 is already close 
to the corresponding result expected in the infinite temperature limit. 
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Figure 1: Temperature dependence of the second order expansion coefficients 
c 2 and c 4 of the pressure in powers of fi q /T obtained in 2-flavour QCD [3]. 
In Fig. la we show c 2 (T) and Fig. lb shows the ratio c±{T) / C2(T) . The 
temperature scale is given in units of the transition temperature at fi q — 
which for the quark masses used in the calculation of Ref . [7j is T ~ 200 MeV. 
For T > To the dashed horizontal line shows the massless Fermi gas value 
of C4/C2. The dashed-dotted curves in (a) show results of a resonance gas 
model calculation for A = 0.9, 1.0, 1.1, 1.2 (from top to bottom) as discussed 
in section 3. The dashed-dotted curve in (b) represents the resonance gas 
value c 4 /c 2 = 0.75. 

2 Also at 0(g 2 ) the high temperature expansion terminates at 0((^-) 4 ). This, however, 
changes in the resumed 0(g 3 ) contribution. The complete expansion up to 0(g 6 lng) has 
recently been presented in 



4 



3 Resonance Gas and Boltzmann Approxi- 
mation 



The analysis of experimental data on the production cross sections of various 
hadrons in heavy ion collisions shows astonishingly good agreement with 
corresponding thermal abundances calculated in a hadronic resonance gas 
model at appropriately chosen temperature and chemical potential j^j. Our 
recent analysis of the equation of state calculated on the lattice at \i q = 
also has shown that a gas of non-interacting resonances can provide a good 
description of the low temperature phase of QCD jHj. We want to extent 
here our analysis to the case of non- vanishing chemical potential. 

The partition function of a resonance gas can be specified through the 
mass spectra for the mesonic and baryonic sectors of QCD, respectively. In 
a non-interacting resonance gas the partition function reads, 



\nZ(T,ii B ,V)= £ lnZ*(T,V)+ £ \nZ^(T,^ B ,V) , (4) 

i £ mesons i g baryons 

where Z^. (Z^J denote single particle partition functions for bosons and 
fermions with mass rrii and fi B — 3/z g is the baryon chemical potential. Here 
the fermion partition function contains the contribution from a particle and 
its anti-particle. The total pressure of the resonance gas is then obtained 
from Eq. ^ and builds up as a sum of contributions from particles of mass 
mj. The dependence of the pressure on the chemical potential at a fixed 
temperature is thus entirely due to the baryonic sector. The contribution, 
p m , of baryons of mass m to the total pressure is given by 



Pm _d_ f°°j 

T ~ 2tt 2 Jo 



/ dkk 2 \n (1 + zexp{-e(k)/T})(l + z^expi-e^/T}) (5) 
Jo 1 J 



where z = exp{/x B /T} is the baryonic fugacity with p, B = 3/x g , d is the 
spin-isospin degeneracy factor and e(k) = \/k 2 + m 2 is the relativistic sin- 
gle particle energy. The pressure may be expressed in terms of a fugacity 
expansion as 



p m d /m N 2 x 



£ {y) E(-i) m r2 K 2 (e m /T) ooBh(wr) 



where K% is the Bessel function. 



(6) 
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In the hadronic phase of QCD the relevant temperatures will always be 
smaller than the transition temperature to the plasma phase determined 
in lattice calculations at \ib = 0, i.e. we are interested in the regime 
T<200 MeV [7j. As the mass of the lightest baryon (mjv) is about five 
times larger than this value, the Bessel functions appearing in Eq. |B] can al- 
ways be approximated by the asymptotic form valid for large arguments, i. e. 
K 2 (x) ~ \Jt\: fix exp(— x). This shows that higher order terms in Eq. IHlare 
suppressed by factors exp(— £(m — hb)/T). As long as {m^ — yU#) > T the 
contribution of baryons to the resonance gas partition function is thus well 
approximated by the leading term in Eq. which constitutes the Boltzmann 
approximation. In this case each baryon/anti-baryon pair contributes to the 
pressure with 

^ = ^(?) 2K2(m/T)cOSh( ^ /T) • (7) 
Thus, the total baryonic contribution to the pressure of a resonance gas reads 

p = F(T) cosh^B/T) , (8) 
where F(T) is defined by 

m^Ys^^K^/T) , (9) 

and the sum is taken over all known baryons and their resonances. 

Current lattice studies of the hadronic phase of QCD with non-vanishing 
chemical potential concentrate on a temperature regime 0.8 T < T < T 
with quark chemical potentials fi q < T. The baryon chemical potential 
Hb — 3/i g thus stays considerably smaller than the nucleon mass. Under 
these conditions the Boltzmann approximation will be applicable and the 
dependence of the pressure on the baryon chemical potential will just be me- 
diated through a multiplicative factor as given in Eq. |S1 We stress that this 
simple relation is independent of details of the mass spectrum as long as all 
fermions are sufficiently heavy. The factorization of the part depending on 
the mass spectrum and that depending on the chemical potential, however, 
also relies on the assumption that interactions among the hadrons and reso- 
nances are negligible. The validity of this assumption can be verified through 
a direct comparison with lattice calculations for fi q > and for T < T . 
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Following Eqs. we can specify the results for the change in pressure, 
the quark number density and quark number susceptibility. In order to com- 
pare the predictions of the resonance gas model with lattice results one still 
needs to perform the Taylor expansion up to the same order as given in Eq. El 
In the Boltzmann approximation we have, 

t = F(T)\cosh^-l]^F(T)(cJ^)\cJ^Y) , 



r F(T)[cosh^-l]^F(T)lc 2 ^j +c~4 Vr 
£ = 3F(T)smh^.F(T)( 2 c 2 (^) + 4c 4 (^) 3 ) 



r , 9F(T) cosh ^ ~ F(T) (^25 2 + 12£ 4 (£2 j j , (10) 

with c 2 = 9/2 and c 4 = 27/8. In the resonance gas model the expansion coef- 
ficients introduced in Eq. El are given by c 2n = c 2 „F(T). We note that ratios 
of these quantities indeed are independent of the resonance mass spectrum 
and only depend on the chemical potential. For fixed \x q jT we thus expect to 
find that any ratio of two of the above quantities is temperature independent 
in the hadronic phase. Using the same order of the Taylor expansion as used 
in the lattice calculations such ratios only depend on c 4 /c 2 = 3/4, i.e. the 
resonance gas model yields a temperature independent ratio c 4 /c 2 . As can 
be seen in the right-hand part of Fig. ^ this is indeed in good agreement 
with the lattice results. We note that this result is independent of details of 
the hadron mass spectrum. It thus should also be insensitive to the change 
in the quark mass used in the lattice calculation. In Fig. |2] we show the 
ratio Ap/T 2 x q for two values of the chemical potential. The good agreement 
between lattice calculations and the hadronic gas results merely reflects the 
agreement found already for the ratio c 4 /c 2 . In the resonance gas model we 
can, however, provide also the complete result for this ratio, 

Al> 1 <\ - cosh-^S/ijT)) . (11) 



T 2 Xq 9 

This is shown as a dashed-dotted line in Fig. El 

The agreement of the ratio c 4 /c 2 calculated in the resonance gas model 
and on the lattice also implies that, for all values of fi q /T, the temperature 
dependence of Ap and its derivatives like e.g. the quark number suscepti- 
bility is to a large extent controlled by the same function, c 2 (T) = \F(T). 
To determine, however, the function F(T) we have to specify the baryon 
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A(P/T 4 )/(X/T S )l u/T=0Bt04 
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Figure 2: The ratio of pressure and quark number susceptibility versus 
temperature for fixed values of fi q /T. The horizontal lines are the results of 
hadron resonance gas model. The points are the lattice values from Ref . [I] . 
While the dashed-dotted curves represent the complete expression (Eq. fTT| 
the dashed curves give the result of a Taylor expansion to the same order as 
that used in the lattice calculations. 



mass spectrum explicitly. This is, in general, known experimentally. How- 
ever, to facilitate a direct comparison with lattice calculations we have to 
take into account that the spectrum is distorted due to the still quite large 
quark masses used in calculations of the equation of state 3 . In addition we 
have to take into account that also the transition temperature is quark mass 
dependent. For the quark mass value used in the numerical study of the 
equation of state jlj [7j the transition temperature has been determined as 
T ~ 200 MeV. While an extrapolation to the chiral limit yields the transi- 
tion temperature of about 170 MeV. We use T to express the hadron masses 
in units of the temperature, m/T = (m/T Q ) ■ (T /T) with T scaled in units 
of T . 

The distortion of the hadron mass spectrum due to unphysically large 

3 The phase transition temperature has been calculated at a large set of quark mass 
values, including rather small values which lead to "almost" physical hadron masses. So 
far the equation of state, however, has been studied with improved actions only for one 
set of quark masses corresponding to a pion mass of about 770 MeV [7]. 
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quark masses (m q ) can in general be deduced from lattice calculations at 
zero temperature. A generic feature of such studies is that the deviation from 
the physical mass value due to unphysically large values of the quark mass 
becomes smaller for heavier hadronic states (see eg. JO]). Moreover, one 
finds [TTl IT2*] that the quark mass dependence is well parametrized through 
the relation, (m#a) 2 = {jn H a) 2 hys + b(m n a) 2 , where (mHCi) p hys denotes the 
physical mass value of a hadron expressed in lattice units and (mjfl) is 
the value calculated on the lattice for a certain value of the quark mass or 
equivalently a certain value of the pion mass (m 2 ~ m q ) . Until now, however, 
the masses of only a few baryonic states constructed from (u, <i)-quarks have 
been studied in more detail on the lattice [TOJ HU H2] • This is obviously not 
sufficient to fix the function F(T) in Eq.(J7J) that requires the contributions 
from a large set of baryonic resonances. 

The above quadratic parametrization of the quark mass dependence of 
baryon masses shows at least for nucleon, delta and their parity partners 
only a weak dependence on the hadron mass. We thus take this as a general 
ansatz for the parametrization of the dependence of baryon masses on the 
pion mass, 



m(mj] . m 2 



l + A^- , (12) 
inn ™>h 

where m(m n ) is the distorted hadron mass at fixed m n and m# is its corre- 
sponding physical value. This parametrization is consistent with our previous 
analysis 6\ where we have used the MIT bag model to determine the in- 
dependence of hadron masses. 

The lattice results [3] for QCD thermodynamics at finite \i q were ob- 
tained in 2-flavour QCD. An immediate consequence of the restriction to 
only two quark flavour is that we have to suppress the contribution from 
strange baryons to the resonance gas in the low temperature phase. More- 
over, since the lattice results were obtained with a quite large value of m q , 
corresponding to m n / ' yfa = 1.84 ± 0.04, we account for modifications of the 
baryon mass spectrum using Eq. (|12jl. From the pion mass dependence of the 
nucleon, delta and their parity partners E] we estimate 0.9 < A < 1.2 
in Eq. ()12J) . This range of values is also expected from the bag model study 
inRef. 0. 

The above discussion fixes our parametrization of the baryonic sector of 
the resonance gas model at unphysical values of the quark mass as they are 
used in current lattice calculations. The resulting temperature dependence 
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of the expansion coefficient Oi(T) is shown in Fig. ^i. The corresponding 
result for the quark number susceptibility at different values of the quark 
chemical potential is shown in Fig. El for the choice A = 1.0. The agreement 
of the resonance gas model and results obtained from lattice calculations is 
indeed quite satisfactory. This indicates that the thermodynamics of the 
confined phase of QCD at finite density is to large extent governed by the 
baryonic resonances. As can be seen in Fig. the 15% uncertainty in the 
parametrization of the baryonic mass spectrum (Eq. IT2jl results in 20% error 
on the values of physical observables, i.e. C2(T), at T = T . 




t/t 



Figure 3: Lattice results from Ref. [3] for the quark number susceptibility 
in 2-flavour QCD calculated in next-to-leading order Taylor expansion for 
different values of the quark chemical potential. The lines are results from 
the resonance gas model using a distorted baryon spectrum (Eq. ED with 
A = 1) and treated within the same approximation as in the lattice study. 



4 Relaxing the lattice constraints 

We want to discuss here in somewhat more detail what the resonance gas 
model calculation suggests for the quark mass dependence of current thermo- 
dynamic studies on the lattice and what the influence of the truncation of the 
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Taylor series expansion on the behavior of thermodynamic quantities in the 
hadronic phase could be. The latter clearly depends on the observable under 
consideration. As anticipated also in Ref. jl] the influence of a truncation of 
the Taylor expansion for the pressure at 0((fi q /T) A ) becomes more severe in 
calculations of the quark number susceptibility as the expansion stops here 
already at 0{{fi q /T) 2 ). The resonance gas model suggests that for \i q jT = 1 
the truncated result for the pressure differs only by 15% from the complete 
result. These truncation errors rise to about 80% in calculations of Xq/T 2 
at fiq/T = 1. The major part of this truncation error could be removed 
by calculating the 0((fi q /T) e ) contribution to Ap/T 4 . These properties are 
seen in Fig. 0] where the results of the resonance gas calculation performed 
with the complete (/i q /T)-dependence and a Taylor expansion truncated at 
0((/i g /T) 4 ) are shown for pressure and quark number susceptibility. 

Already in the discussion of the Ap/T 2 x q we indicated that within the 
Boltzmann approximation this and similar ratios are independent of the res- 
onance mass spectrum and thus also on the quark masses. Changes in the 
quark mass thus will influence calculations of the pressure and its derivatives 
in a similar way. Replacing in the resonance gas calculation the modified 
baryon spectrum by the experimentally known spectrum will increase the 
value of the pressure as all baryons become lighter. This effect is somewhat 
reduced as also the relevant temperature scale is shifted to smaller values, i. e. 
the transition temperature will shift from T ~ 200 MeV to T ~ 170 MeV. 
As can be seen in Fig. 0] this increases the pressure and its derivatives. The 
extrapolation to the physical case depends, however, quite sensitively on the 
value for the critical temperature. We note that the resonance gas model 
calculation favours a small critical temperature as it seems to be unlikely 
that Ap/T 4 will exceed the corresponding ideal gas value at T . 

Finally, we want to comment on the convergence radius of the Taylor 
expansion in (fi q /T). The resonance gas model, of course, does not lead to 
critical behaviour. Consequently also the dependence on (/i g /T) is given by 
an analytic function. The resulting Taylor expansion has an infinite conver- 
gence radius, which in terms of the convergence criterium used in Ref. 0] is 
reflected in the fact that, 



for all temperatures in the resonance gas model. In the case of QCD, we 
expect, however, that the convergence radius is bounded from above by the 
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Figure 4: Quark number susceptibility (a) and change in pressure (b) for 
fixed quark chemical potential fi q /T = 1 as a function of T/To. The points 
are lattice results from Ref. 0] and lines are the resonance gas model re- 
sults. The dashed and dashed dotted lines are obtained with a baryon mass 
spectrum appropriate for the unphysical quark masses used in lattice calcu- 
lations and with a Taylor expansion truncated at C((/i g /T) 4 ) and the full 
result, respectively. The full lines are resonance gas model results obtained 
with physical hadron masses, no expansion in fi g /T and for three values of 
the transition temperature, T = 160 MeV (lower), 170 MeV (middle) and 
180 MeV (upper), which cover the range of current lattice estimates for the 
chiral limit extrapolation of T c in 2- flavour QCD. 

location of the phase boundary to the quark gluon plasma phase. In partic- 
ular, we expect that the ratios c 2n /c 2n+2 stay close to unity for temperatures 
close to Tq. In fact, this also is the case for the low order expansion coeffi- 
cients in the resonance gas model, i.e. c 2 /c4 = 4/3. This, however, changes 
already at the next order, c^/cq = 10/3. We thus expect that differences be- 
tween lattice results on the QCD thermodynamics and resonance gas model 
calculations should show up at 0((fi q /T) 6 ). 
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5 Conclusion and Outlook 



We have shown that basic features of the density dependence of the QCD 
equation of state in the hadronic phase observed in recent lattice studies can 
be understood in terms of the thermodynamics of a baryonic resonance gas. 
A quite robust result, independent of the detailed structure of the hadron 
mass spectrum, is the dependence of thermodynamic quantities on [i q jT at 
fixed temperature. This indicates that the thermodynamics at low temper- 
atures is dominated by heavy degrees of freedom which justify a Boltzmann 
approximation for the partition function. Within the resonance gas model 
these degrees of freedom are non-interacting which leads to simple relations 
among different thermodynamic observables. In particular, we argued that 
the ratio Ap/T 2 x q is independent of temperature at fixed \i q jT . This is in 
agreement with current lattice calculations. It is, however, clear that this 
has to change, when the system undergoes a true phase transition at some 
temperature T c (fi q ) < T for a sufficiently large value of \i q jT . In particu- 
lar, we expect that the quark number susceptibility will diverge at the chiral 
critical point, i.e. at the second order endpoint of a line of first order phase 
transitions [13J, which should lead to a dip in the ratio Ap/T 2 x q at T c (/i q ). 

Finally we note that the relevance of the resonance gas model predictions 
for the description of the low temperature phase of QCD can also be tested 
in the imaginary chemical potential approach ^3] which also is used in lat- 
tice calculations. The resonance gas model, for instance, suggests that the 
change in pressure is quite well described by a simple trigonometric function, 
cos(3// g /T). 
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